Kondo effect in Normal-Superconductor Quantum Dots 
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Abstract 



We study the transport properties of a quantum dot coupled to a normal and 



a superconducting lead. The dot is represented by a generalized Anderson 
model. Correlation effects are taken into account by an appropriate self- 
energy which interpolates between the limits of weak and strong coupling 
^SJ ■ to the leads. The transport properties of the system are controlled by the 

interplay between the Kondo effect and Andreev reflection processes. We 
show that, depending on the parameters range the conductance can either 
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be enhanced or suppressed as compared to the normal case. In particular, 
by adequately tunning the coupling to the leads one can reach the maximum 



value 4e 2 //i for the conductance. 
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The Kondo effect is a prototypical correlation effect in Solid State Physics. Although it 
was first analyzed for the case of magnetic impurities in metals, in the last years there has 
been a renewed interest in Kondo physics with its observation in a semiconductor quantum 
dot (QD) |],@]. Quantum dots constitute an ideal laboratory for testing the theoretical 
predictions as they allow to vary the relevant parameters in the problem in a controlled way. 
This technology also opens some new possibilities like the exploration of the Kondo effect 
when the dot is connected to superconducting leads. The interesting issue in this case is 
related to the competition between the strong Coulomb interaction in the quantum dot and 
the pairing interaction within the leads. 

From the theoretical side, the Kondo effect in QDs has been mainly analyzed by means 
of the single- level Anderson model H^. The theory predicts an enhancement of the dot 
conductance at low temperatures due to the development of the so-called Kondo resonance. 
The case when one of the leads is superconducting has been recently analyzed by some 
authors using different theoretical methods 0-0] assuming a modified Anderson model in 
which one of the metallic electrodes is substituted by a BCS superconductor. While some 
authors have predicted an enhancement of the conductance due to Andreev reflection at 
the superconducting lead ||, others have predicted the opposite effect j7|. In Refs. |5],[7| the 
infinite charging energy limit (U — » oo) has been assumed. However, in an actual experiment 
this assumption may not be completely justified (for instance, in the experiments of Ref. 
U the ratio U/T, T being the dot tunneling rate, was estimated to be around 6.5). The 
approach presented in this letter would allow to analyze this problem for a broad range of 
the different parameters of the model. We will show that, depending on the values of these 
parameters, one can obtain either an enhancement or a reduction of the conductance with 
respect to the normal case. 

Our approximation scheme is based on the hypothesis that a good approximation to the 
electron self-energy can be found by interpolating between the limits of weak and strong 
coupling to the leads. This interpolative method has been applied successfully to analyze 
different strongly correlated electron systems like the equilibrium |§||, the non-equilibrium 



and the multilevel jL(J Anderson models and the Hubbard model fLl,12|. In this letter 



we shall discuss how to extend this method to the superconducting case. 
For describing a N-QD-S system we use an Anderson-like Hamiltonian 

H = H N + H S + Y, ^K + Un ] h i + H T , (1) 

(7 

where n a = d\d a , H N and H s represent the uncoupled normal and superconducting leads 
respectively; Ht = J2keN Sa^kd^c^a + h.c. describing the coupling between the dot level 
and the leads. Within this model the dot is represented by a single spin degenerate level 
with a repulsive Coulomb interaction described by the U-term in Eq. (1). We shall assume 
that the superconducting lead is well described by the BCS theory with a superconducting 
gap A and the normal lead is, as usual, characterized by a flat density of states around the 
Fermi level, pf- 

The transport properties of this model can be obtained by means of Green function 
techniques. In order to analyze the linear regime the main quantity to be determined is the 
dot retarded Green function, which in a Nambu 2x2 representation adopts the form 



G r (u) 



Lul-e a z - S r '(u;) -T N (u) -T s (uj) , (2) 
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where f N and f s are the tunneling rates given by t N = T L I, and the superconducting 
tunneling rate f 5 is given by f s = T R g, where T L , R = id\ R p F , 9\\ = 9ti = -uj/y/A 2 - u 2 



and <7i2 = #21 = A/VA 2 — u 2 (the chemical potential of the superconducting lead is taken 
as zero). The self-energy £ r (cu) takes into account the effect of Coulomb interactions. To 
the lowest order in U this is given by the Hartree-Fock Bogoliubov approximation: S r = 
U < n > & z + A^iXr, Arf being the proximity effect induced order parameter in the QD, 
A<2 = U(dld\). The crucial problem is to find a good approximation to include correlation 
effects beyond this mean field approximation. 

Within the spirit of the interpolative method commented above, the self-energy is con- 
structed in such a way as to interpolate between the limits of weak and strong coupling to 
the leads, for which the exact result are known. Let us first analyze the weak coupling or 



atomic limit. In this case we have ttfi/U —* and thus f n/U, Ts/U —>■ 0. In this limit the 
induced order parameter in the QD vanishes faster than (tji/U) 2 and one can neglect the 
non-diagonal elements in the self-energy matrix. On the other hand, the diagonal elements 
can be easily evaluated in this limit using the equation of motion method || and have the 
form 

E; i22 ^ ± C/<n> + t/2<ft ^- <fl> > (3) 

u ' iZ w^e =F U(l- < n >) 

In the opposite limit, U/t LR — > 0, one can accurately evaluate the self-energy using stan- 
dard perturbation theory in the Coulomb interaction. The different diagrams contributing 
to the second order self-energy are depicted in Fig 1. In the superconducting case, there 
appear additional diagrams to the one in the normal case (diagram a) corresponding to the 
interaction of an electron with an electron-hole pair in the QD; the remaining diagrams 
contain at least one anomalous propagator and vanish identically in the normal state. As in 
the normal case ||, the non-perturbed one-electron Hamiltonian, over which the diagram- 
matic series is constructed, is taken as an effective mean field, characterized by an effective 
dot level e e //, having the same dot charge as the fully interacting problem. As shown in 
Ref. H this self-consistency condition provides in the normal case a good fulfillment of the 
Friedel sum rule at zero temperature. The extension of this procedure to the superconduct- 
ing case requires dressing the propagators in the diagrams of Fig. 1 with the non-diagonal 
self-energy £12 in order to impose also consistency in the non-diagonal charge (d^d^). Notice 
that although the interaction in the QD is repulsive, there is always some induced paring 
potential in the dot due to the proximity effect. The inclusion of this effect for finite U is 
very important for the correct description of the dot electronic properties. 

The original interpolative scheme stems from the observation that the second order self- 
energy (S^ 2 )) has a similar functional form as the atomic self-energy for large frequencies 
H thus allowing for a smooth interpolation between the two limits. In the superconducting 

case, the diagonal elements of the second order self-energy behave as 

(2) U 2 <n> (l-<n>) 

^11,22 ~ — l 4 J 

w =F e eff 
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for large frequencies, while the non-diagonal elements decay faster than U 2 Juj. This behavior 
permits to define a Nambu 2x2 interpolative ansatz for the self-energy matrix as: 

E(w) = U(n)a z + A d a x + [i - at {2) a z ] ' ± {2) (cu), (5) 

where 

e + (1- < n >)U -e e tf 

a = — . 

U 2 <h> (1- < n >) 

and E*- 2 -* is the second order self-energy matrix whose elements are given by the diagrams 
depicted in Fig. 1. 

Using this ansatz one recovers the correct behavior of the self-energy both in the weak 
and strong coupling limits. Moreover, this ansatz satisfies the exact relations between the 
different matrix elements, i.e E 12 (^) = S 21 (c<j) and En (a;) = — Ej^— ui). 

Due to the presence of an additional energy scale fixed by A the number of different 
physical regimes is larger than in the normal case. We will mainly consider the more inter- 
esting physical regime r = Tl + Yr ~ A |TJ| . In Fig. 2a we show the dot spectral density 



(LDOS) for a symmetric case (eo = —U/2) with Tl = Tr ■ = A and increasing values of 
U. As can be observed, when U < A the LDOS exhibits a double peak around the Fermi 
energy which is due to the influence of the superconducting electrode by the proximity effect. 
However, as U increases the double peak is replaced by a single narrow Kondo resonance as 
in the normal case. The comparison with the normal case reveals that the Kondo resonance 
gets narrower in the superconducting case and its height increases with U above the normal 
value. In the limit U — > oo, this height approaches the value 2/(irT), which is twice the value 
in the normal case at zero temperature, as fixed by the Friedel sum rule. The narrowing of 
the Kondo resonance gives rise to a lowering of the Kondo temperature with respect to the 
normal case. 

For energies larger than A the differences between the normal and the superconducting 
LDOS become negligible, with the usual broad resonances at e and e + U which become 
more pronounced for increasing U. 



By varying the dot level position eo one can study the transition from the Kondo to 
the mixed valence regime. The evolution of the dot LDOS is illustrated in Fig. 2b. When 
approaching the mixed valence regime (|eo| < T or |eo + U\ < Y) the Kondo resonance is 
replaced by an asymmetric broad resonance close to the Fermi energy as in the normal case. 
In the superconducting case, however, the LDOS develops an additional structure associated 
with the BCS divergencies at the gap edges @. 

As in any NS contact, transport at low voltages is possible due to Andreev reflection 



processes. At finite temperature, the linear conductance is given by the expression |L4 



16e 2 r°° ( df\ 

G = —Y L J JE\m (G r 12 G a u ) (Y R - ReS 12 ) l-^L) 



(6) 



where f(E) is the Fermi function. At zero temperature, ImE(O) = 0, and Eq. (^ reduces 

to 

4e 2 4r?f 2 
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where f R = Y R — ReS 12 (0) and e = e + ReS n (0). Notice that Eq. (|7|) coincides at U = 
with the well known non-interacting result | |15| . 

One would expect that for a dot symmetrically coupled to the leads (i.e. Tl = T R ) and in 
the case of electron-hole symmetry (eo = —U/2), the conductance should reach its maximum 
value 4e 2 /h 0. However, the actual situation is more complex due to the reduction of the 
induced paring amplitude in the dot arising from the repulsive Coulomb interaction. As a 
consequence the conductance decreases for increasing U even in this case. This decrease is 
illustrated in Fig. 3a where we plot the conductance as a function of U in the symmetric case 
for different values of T/A. For large U/Y we find that the conductance decreases roughly 
as (Y/U) A . This behavior can be understood as follows: in order to have a vanishing pairing 
amplitude in the U/Y —>■ oo limit, the non-diagonal self-energy S 12 should tend to cancel 
the non-diagonal tunneling rate (r,s)i2- By analyzing the expression of diagram d in Fig. 1, 
this requires that Gu decays as (Y/U) 2 and therefore the conductance given by Eq. (g) in 
our approximation should decay roughly as (Y/U) 4 . This decay is probably less pronounced 
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than in the exact solution where one would rather expect an exponential behavior in the 
Kondo regime. 

Although the previous analysis shows that the maximum value for the conductance 4e 2 /h 
can never be reached in the symmetric case for finite U, this is not necessarily the case for 
an asymmetric situation with T L ^ Y R . In fact, if the coupling to the electrodes could be 
tunned in order to reach the condition Tr = Tl then, Eq. @ predicts a maximum in the 
value of G. As shown in Fig. 3b, this condition can be reached by increasing the coupling 
to the superconducting electrode. The ratio between Tr and Tl at the maximum becomes 
larger for increasing U . In a situation with electron hole-symmetry, like the one depicted in 
Fig. 3b, the conductance at zero temperature reaches its maximum possible value 4e 2 /h. 

In normal quantum dots a signature of the Kondo effect is given by an anomalous tem- 
perature dependence in the linear conductance |IJ, which exhibits a continuous transition 
from a maximum conductance in the Kondo regime to well resolved conductance peaks as- 
sociated with Coulomb blockade. When one of the electrodes is superconducting there is 
also a decrease of conductance with temperature in the Kondo regime. However, as depicted 
in Fig. 4, the conductance already exhibits a double peaked structure at zero temperature 
when T^ = T/j. The reduction of conductance with temperature is in this case much faster 
than in the normal case, as shown in Fig. 4 (inset). This difference is a consequence of the 
lowering of the Kondo temperature due to the presence of the superconducting electrode. 

In conclusion, we have analyzed the electronic transport properties of a quantum dot 
coupled to a normal and a superconducting lead. For this purpose we have introduced an 
electron self-energy which interpolates between the limits of weak and strong coupling to 
the leads, an approach which has been previously used for normal systems f4]j8|-|12|| . This 
approximation allows to describe a broad range of parameters including the relevant one for 
an actual experiment. On the other hand, we have shown that for finite charging energy 
the dot conductance can either be enhanced or suppressed with respect to the normal case. 
While in a symmetrically coupled dot (r^ = T R ) an increasing charging energy tends to 
reduce the conductance, in the asymmetric case it is always possible to reach a maximum in 



the conductance by fine tunning the coupling to the superconducting electrode. In the case 
of electron-hole symmetry this maximum reaches the value 4e 2 //i at zero temperature. The 
predictions presented in this work could be tested experimentally using similar technologies 
to those currently used for normal quantum dots fI|,[^,|T6|1 
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FIG. 1. Second order self-energy diagrams. 
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FIG. 2. a) Dot spectral density in the symmetric case (eo = —U/2 and Tl = Tr) with T = A 
for different values of U/T. The inset shows a blow up of the region around the Fermi energy, b) 
Dot spectral density for different values of eo with U = 10T and T = A. 
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FIG. 3. a) Conductance at zero temperature for the symmetric case as a function of U/T and 
for different values of T/A. From bottom to top T/A = 0.125, 0.25, 0.5, 1.0, 2.0, 4.0, 8.0. b) same 
as (a) for asymmetric coupling to the leads as a function of Tr/Tl and different values of U/Tl- 
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FIG. 4. Conductance for different temperature values: U = WT, A = T/2 and T/cf = 0.0, 
0.0005, 0.001, 0.0025, 0.005, 0.01. Inset: normalized conductance as a function of temperature for 
N-dot-S (full line) and N-dot-N (dashed line) at E F = U/2. 
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